The Deformation of an Elastic Substrate by a Three-Phase Contact Line 
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Young's classic analysis of the equilibrium of a three-phase contact line ignores the out-of-plane 
component of the liquid-vapor surface tension. While it has long been appreciated that this un- 
resolved force must be balanced by elastic deformation of the solid substrate, a definitive analysis 
has remained elusive because conventional idealizations of the substrate imply a divergence of stress 
at the contact line. While a number of theories of have been presented to cut off the divergence, 
none of them have provided reasonable agreement with experimental data. We measure surface 
and bulk deformation of a thin elastic film near a three-phase contact line using fluorescence con- 
focal microscopy. The out-of-plane deformation is well fit by a linear elastic theory incorporating 
an out-of-plane restoring force due to the surface tension of the gel. This theory predicts that 
the deformation profile near the contact line is scale-free and independent of the substrate elastic 
modulus. 



At first glance, it is hard to imagine that there is any- 
thing mysterious about a droplet of water resting on a 
solid surface. However, mathematical descriptions of an 
idealized contact line, where the liquid, solid and vapor 
phases meet, can contain perplexing singularities. For 
example, the diffusion-limited evaporation rate of a sta- 
tionary droplet [l| and the strain-rate of a translating 
droplet [3, y[ both diverge at the contact line. 

As articulated by Young, the equilibrium angle of a 
three-phase contact line is determined by a balance of 
solid-vapor, solid-liquid, and liquid-vapor surface ten- 
sions 4]. Surface tension is the derivative of the free 
energy with respect to area and can be visualized as a 
generalized force per unit length acting on the contact 
line, as shown in Fig. [TJ Specifically, the solid-liquid 
surface tension, 7 S ;, exerts a force per unit length ra- 
dially inward, and the solid- vapor surface tension, j sv , 
exerts a force per unit length radially outward. The con- 
tact angle sets itself so that the horizontal component 
of the liquid- vapor surface tension, cos(0)7z„, balances 
7s; — Jsv However, this classic force-balance leaves the 
vertical component of the liquid-vapor surface tension 
unbalanced. While it is clear that the normal force due 
to the liquid-vapor surface tension must be balanced by 
elastic deformation of the solid substrate, the calculation 
of the deformation poses important conceptual difficul- 
ties. In particular, the force per unit area, or stress, ex- 
erted by the liquid- vapor surface tension diverges at the 
contact line. Continuum elastic theory therefore predicts 
a strain divergence at the contact line. 

To avoid this difficulty, Shanahan and DeGennes im- 
posed a length-scale cut-off and focused their analysis on 
regions far from the contact line 0, Alternatively, 
Long et. al. 7\ suggested that the divergence in the 
strain could be cut off by including a surface-tension 



vapor 




FIG. 1: Schematic diagram showing a sessile droplet on an 
elastic substrate. Here, 7 S ; = j av . Superimposed on the right 
side of the droplet is the conventional force balance used to 
derive Young's equation. The anticipated scale-free deforma- 
tion of the substrate at the contact line is superimposed on 
the left side. 



penalty for the additional surface area of the deformed 
solid. Pericet-Camara et. al. Q recently measured the 
topography at the free surface of a PDMS substrate due 
to a drop of ionic liquid. Their data showed agreement 
with Long's theory in the long and short wavevector lim- 
its, but large discrepancies were found at distances from 
the contact line comparable to the substrate thickness. 

In this Letter, we report measurements of the in-plane 
and normal displacement fields in a gel near its contact 
line with a sessile droplet of water. We find that a linear 
elastic theory accounting for the finite thickness and sur- 
face tension of the gel resolves some of the conceptual dif- 
ficulties imposed by the stress singularity and provides an 
excellent quantitative description of the out-of-plane de- 
formation of the gel over all accessible length scales. In- 
terestingly, the deformation near the contact line is scale 
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FIG. 2: Deformation of a thin elastic substrate due to a three-phase contact line. Blue points indicate deformations u z (a;) (top) 
and «! (i) (bottom) measured at the surface (z = h, left) and just below the surface (h — z = 3 /im, right) of a gel near a 
three phase contact line. The vapor phase is to the left of the peak and the liquid drop lies to the right of the peak. Best fits 
to various theoretical models are superimposed on the data. The dashed lines are for a finite thickness elastic slab and the 
solid lines are for a finite thickness elastic slab with surface tension. a;-shaped plot symbols for u x (h) indicate points subject 
to systematic errors due to removal and re-deposition of beads by the contact line. 



free and depends only on the ratio of the solid and liq- 
uid surface tensions. While a careful accounting for sur- 
face tension gives a good description of the out-of-plane 
displacements, discrepancies remain for the in-plane dis- 
placements. 

We measure the deformation of a highly-elastic sili- 
cone gel (Dow Corning Toray, CY52-276A/B) due to a 
10 //L droplet of de-ionized water. Uncured gel was spin- 
coated onto a relatively rigid glass coverslip to form a 
film of thickness ft, = 20 /im. We estimated the Young's 
modulus of the gel to be 3 kPa using bulk rheology. To 
quantify the deformation field, two layers of fluorescent 
beads were incorporated into the film, one at the free 
surface and one embedded 3 fj,m below the surface @- 
A spinning disc confocal microscope (Andor Revolution, 
mounted on a Nikon Ti Eclipse inverted microscope with 
an oil- immersion 40X objective, NA = 1.3) was used to 
image the beads. Three-dimensional image stacks of the 
fluorescent tracer particles were acquired for a 5 minute 
period beginning about 15 minutes after the droplet was 
deposited on the surface. Over this period, the contact 
line moved smoothly at approximately 36 /tm/min as 
the drop evaporated. An image stack was also acquired 
about 45 minutes later to provide a zero-stress reference 
for bead positions. Bead displacements were determined 
using centroid analysis and particle tracking software in 
MATLAB [l(J El ■ Since we expect forces to be invari- 
ant along the contact line, we average the bead displace- 
ments in this direction. On the assumption that the force 
exerted by the contact line did not change significantly 
from one timestep to the next, we combined lab-frame 



displacement profiles from four successive timesteps sep- 
arated by 55 sec to construct a displacement profile as 
a function of distance from the contact line, as shown in 
shown in Fig. [2] 

The solid substrate forms a symmetric ridge about 1 
/im high just below the contact line. A wide and shallow 
valley appears on either side of the ridge, with resolvable 
deformations observable up to 60 /im from the contact 
line. In-plane displacements point toward the contact 
line and decay very slowly. There is a slight but un- 
mistakeable asymmetry in the ^-displacements on either 
side of the contact line. This is surprising because the 
residual force of a static contact line on a smooth surface 
should be entirely out of plane, even if the contact angle 
is not identically 90°. We suspect that this small asym- 
metry is due to either pinning forces or viscous stresses 
near the slowly receding contact line. As expected, the 
displacement fields are bounded, with no hint of a strain 
singularity near the contact line. 

To construct a model that explains these observations, 
we begin with linear elastic theory. The linearized gov- 
erning equation for the displacement field in an isotropic 
elastic solid is:(I3| 

(1 - 2v) V 2 m + V(V ■ u) = 0, (1) 

where v is Poisson's ratio (v = i for incompressible ma- 
terials), and any forces are exerted at the boundaries. 
Boussinesq solved this equation in the case of an infinitely 
thick substrate [l2j]. Integrating this solution along a line, 
we find the extensional strain, e zz , at a depth d directly 
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below the contact line in an incompressible material, 
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where 7/„ is the surface tension of the liquid as indicated 
in Fig. [TJ and E is the Young's modulus of the substrate. 
For simplicity, we assume a contact angle 9 = ir/2 and 
that d is much less than the length of the line. Thus, 
strain diverges as e zz ~ 4 approaching the contact line. 

While Boussinesq's solution assumes a semi-infinite 
slab of elastic material, h/R 1, our gel is much thinner 
than the radius of the drop, h/R « 10~ 2 . We recently 
presented a solution to Eq. ([I} for a finite-thickness elas- 
tic film, which accounts for the zero-displacement bound- 
ary condition at the interface with a rigid substrate, lo- 
cated, in our coordinate system, at z = [9j. In Fourier 
space, the stresses at the surface, <Ti Z (h), are linearly re- 
lated to displacements at height z, Uj(z): 
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where the tensor Q is a generalized spring constant and 
we sum over repeated indices. For a three phase con- 
tact line, we exploit translational invariance along the 
contact line and work in two dimensions. Assuming 
that residual force from the contact line is normal to 
the surface, the stress due to the line can be written: 
Ciz(fc) = Sizjiv sin 9. Using the solution presented in Ref. 
[9(, the z-displacements at the surface of an incompress- 
ible material due to the contact line are given by: 
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Where k is the wavenumber in the x-direction. This 
solution reduces to the Boussinesq form in the limit 
of h — y oo: Q zz = ^grCT- As shown in Fig. [3j 
two salient features differentiate the finite thickness so- 
lution (blue dashed line) from the Boussinesq solution 
(red dashed line). First, the substrate thickness h in- 
troduces a length-scale so that Q~ z has a maximum at 
around k — 1/h. This leads to the dimples on either 
side of the peak. Second, the finite thickness solution 
suppresses a divergence in Boussinesq's solution as \k\ 
goes to zero. This ensures that displacements go to zero 
as x goes to ±oo. However, the divergence at the ori- 
gin remains. Mathematically, the local behavior of the 
real-space z— displacement at x — is governed by the 
inverse Fourier transform of lim Q^ z ■ In this limit, the 

\k\ — >oo 

finite thickness solution and the Boussinesq solution are 
equivalent: in both, u z diverges logarithmically, and e zz 
diverges like 1/z. Physically, the thickness is expected to 
be irrelevant in this region, since the solution should be 
dominated by local quantities near the singularity. 

The measured displacement profiles were fit with this 
solution. Since the film thickness, h, and the positions of 
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FIG. 3: Fourier transform of the deformation of an elastic 
substrate due to a out-of-plane line force. The red dashed 
line shows the prediction of the classic theory of Boussinesq, 
which decays like 1/k. The blue dashed line supresses the 
divergence at k = by accounting for the finite thickness of 
the substrate. The black line incorporates the surface tension 
of the substrate and decays as 1/k 2 at large k. 



the layers of beads, z, are known, this leaves one free fit- 
ting parameter, r yi v sin9/E. We found the best-fit using 
least-squares minimization of u z for both layers of beads 
and u x for the inner layer of beads simultaneously. Since 
the contact line removed and redeposited beads on the 
surface, we did not include u x from this layer in the fit. 
The best fit is shown by the blue dashed line in Fig. [5] 
This theory not only fails to predict a bounded deforma- 
tion at the surface, but it also places the minima in u z too 
close the contact line, and systematically underestimates 
the in-plane deformations. 

To improve this theoretical prediction, we take into 
account the surface tension of the gel and introduce a 
free energy penalty for the creation of new surface area on 
the substrate. Following Long Q, we continue to assume 
7s; = Isv = Is, which is valid for an equilibrium contact 
angle 9 = n/2. The measured equilibrium contact angle 
was about 105°, so this approximation is expected to be 
reasonable but not completely accurate. 

Since the surface tension of the solid exerts a force per 
unit length like the surface tension of the liquid, it is 
capable of cutting off the singularity at the origin. While 
Eq. ([3]) still governs displacements within the bulk, the 
stress at the surface now includes an internal contribution 
a-y from the substrate surface tension, as well as the term 
due directly to the contact line, a ex t- 
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To arrive at the form of the surface stress, we first 
calculate the free energy per unit length due to ad- 
ditional surface area, \u 2 z k 2r y s . The associated stress 
is the derivative of this free energy with respect to 
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small changes in the z— displacement of the free surface: 
k 2 "fsU z (h)/2, [7[. Thus, we can write: 

<Jext,iz = QSij(z,h)uj(z), (6) 

where 

QSij(z,h) = Qij(z,h) + -^ — -Qik(h,h)Q k: j{z,K) (7) 



This result, which incorporates the solid surface tension 
into the finite-thickness solution, is plotted as a solid 
black curve in Fig. [3] While it agrees with Eq. (j4]) 
for small wavenumbers, it decays much faster at large 
wavenumbers, QS~^ ~ l/|fc| 2 . Importantly, this fast de- 
cay ensures that u z and e zz remain finite at the contact 
line. Here, the shape of the surface is a cusp given by 
u z (0) — u z (x) — ~fi v sin0|x|/7 s . This self-similar shape 
applies in the regime \k\ ^> Ej^ s : at these small length 
scales, the elastic modulus does not effect the shape of 
the cusp. Finally, the solid surface tension pushes the 
peak to lower k, increasing the characteristic wavelength 
of the displacement profile. 

This solution provides a much closer fit to the data. 
The best-fit curves are shown as solid black lines in 
Fig. [5] The values extracted for the slope of the cusp 
ji v sinO/js = 0.14 and the scale of out-of-plane defor- 
mation jivS'mO/E = 3.2 fim. With the addition of the 
solid surface tension, our theory beautifully captures the 
out-of-plane deformation, including the shape of the peak 
and locations of minima. However, these values do not 
accurately fit the in-plane displacements. This discrep- 
ancy may be due to in-plane forces due to pinning or 
viscous drag near the contact line. 

In conclusion, we measure the displacement field in a 
linear, elastic substrate due to the contact line of a wa- 
ter drop, and show that the z— displacement profiles are 
well-fit by a model that includes both the finite thickness 
of the substrate and the substrate's surface tension. This 
model solves many of the conceptual problems raised by 
applying the Boussinesq solution to the unresolved force 
of an equilibrium contact line. The finite thickness of the 
substrate ensures that displacements go to zero far from 
the contact line. The substrate's surface tension coun- 
ters the stress singularity at the origin, ensuring that 
the z— displacement and strain at the origin do not di- 
verge. Additionally, work is needed to treat cases where 
the equilibrium contact angle is far from 90° , indicating a 



To determine the z— displacements in response to the line 
force, we examine QS~}. Again, for the specific case of 
z = h and v = \'- 
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large discrepancy in surface energies, "f s i and j sv , which 
we have taken to be identical. This work enables a closer 
look at pinning and viscous dissipation near the contact 
line and provides a basis for the extension of traction 
force microscopy. We acknowledge support from Yale 
College, Unilever and NSF DBI-0619674. 
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